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ABSTRACT 

Context. One challenge for measuring the Hubble constant using Classical Cepheids is the calibration of the Leavitt Law or period- 
luminosity relationship. The Baade-Wesselink method for distance determination to Cepheids relies on the ratio of the measured radial 
velocity and pulsation velocity, the so-called projection factor and the ability to measure the stellar angular diameters. 
Aims. We use spherically-symmetric model stellar atmospheres to explore the dependence of the p-factor and angular diameter cor- 
rections as a function of pulsation period. 

Methods. Intensity profiles are computed from a grid of plane-parallel and spherically-symmetric model stellar atmospheres using 
the SAtlas code. Projection factors and angular diameter corrections are determined from these intensity profiles and compared to 
previous results. 

Results. Our predicted geometric period-projection factor relation including previously published state-of-the-art hydrodynamical 
predictions is not with recent observational constraints. We suggest a number of potential resolutions to this discrepancy. The model 
atmosphere geometry also affects predictions for angular diameter corrections used to interpret interferometric observations, suggest- 
ing corrections used in the past underestimated Cepheid angular diameters by 3 - 5%. 

Conclusions. While spherically-symmetric hydrostatic model atmospheres cannot resolve differences between projection factors from 
theory and observations, they do help constrain underlying physics that must be included, including chromospheres and mass loss. 
The models also predict more physically-based limb-darkening corrections for interferometric observations. 

Key words. stars:atmospheres / stars:distances / stars: variable: Cepheids 



1. Introduction 

Classical Cepheids form an important rung on the cosmic 
distance ladder because they follow a period-lum inosity re- 
lation, also called the Leavitt Law jLeavittl Il908h . As such, 

Cepheids are used to measure distances to other salax- 
f i I 1 

ies ( e.g. Pietrzvriski et al. 201 0!), and cosmol ogic al parame- 
ters dFteedman et al. 2001; Riess et al.l |2009l l20Tl) , but the 
Leavi tt Law must be calibrated using Galactic (e.g. iTumerl 
[2010^ and Large Magellanic Cloud (LMC) Cepheids (e.g. 
,Ngeow & Kanbur 2008; Ngeow et al. 2009). 

The Leavitt Law is calibrated by measuring Cepheid dis- 
tances using numerous of techniques, such as parallax, clus- 
ter membership, and the Baade-Wesselink method. The Leavitt 
Law has been derived by assuming all LMC Cepheids are 
at the same distance, where the distance to the LMC is de- 
rived using a nurn b er of methods (e.g. iGieren et aP l2005t 
IClem ent et al.' 2008; 'Bona nos et al.l 1201 ih . Parallax ha s been 
measured for a number of G alactic Cepheids (|van Leeuwe n et alj 
l2007t [Benedict etld] l2007h . but this method is limited by the 
precision of Hipparcos and the Hubble Space Telescope. Only 
a small number of Cepheids have precise parallaxes measured, 
and this will only change when th e GAIA satellite is launched 
(Windmarket al. 2011). Likewise. iTurnerl (2010) measured dis- 
tances to 24 Cepheids belonging to clusters and groups. On 
the other hand, the Baade-Wesselink method can be used for 



Galac tic and LMC Cepheids dKervella et alJl2004l: lOieren et all 

'2005), allowing for a much greater sa mple siz e. 

The Baade-Wesselink method (jBaadd 119261 ; IWesselinM 

Il946h measures the distance to a Cepheid by measuring the 
change of angular diameter and the actual change of radius as 
a function of phase. These measurements yield the mean radius 
and mean angular diameter, hence the distance. The change of 
angular diameter can be determined using a number of meth- 
ods. One method is from interferometric ob servations, where the 
angular diameter can be directly resolved dKervella et al.ll2004l 
2006; Merand et al. 2006). The method uses angular diameter 
corrections for limb darkening when the observations have lim- 
ited resolution. A se cond method is using the infrared surfac e 
brightness technique dBames & Evansll976HGieren et al.ll989l) . 
This method uses a correlation between the angular diameter at 
any phase and the star's V-band flux and color This method is 
derived from the relation FboI °^ ^^^df where 9 is the angu- 
lar diameter, FboI is the bolometric flux and T^ff is the effec- 
tive temperature. The effective temperature is correlated to the 
color based on observations of other stars or model stellar atmo- 
spheres, while the bolometric flux is replaced by the V-band flux 
and the bolometric correction. This method assumes the radius 
is defined at a specific location in the stellar photosphere, i.e. at 
an optical depth r = 2/3. 

The change of radius is AR - J v^uisdl, where Vpuis is the 
pulsation velocity. One cannot directly measure the pulsation ve- 
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locity but instead only the radial velocity. The measured radial 
velocity is the projection of the pulsation velocity weighted by 
the intensity at each point on the stellar disk, i.e. 
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The intensity, Up), is given as a function of the cosine of the an- 
gular distance between a point on the stellar disk and the center 
of the disk, defined as /i. The quantities "TY and K are the sec- 
ond and third moment of intensity, respectively (Mihalas 1978). 
[Getting (Il934 defined the ratio "HIK to be the projection fac- 
tor or p-factor, and it i s a func t ion of the intensity profile of a 
Cepheid at any phase. iGettind (Il934l) determined a projection 
factor p - 1.41 based on an assumed limb-darkening profile. 

This p-factor was the accepted value until Parsons (1972) 
measured it by fitting observed spectral lines using line bisec- 
tors to compare with theoretical spectra. By measuring the line 
center and the doppler shift from the line center of a static 
atmosp here, one can determine the projection factor Parsons 
(119721) measured p = 1.31. Hindsley & Bell (1986) used a 
similar method to measure how the p-factor varies as a func- 
tion of pulsation p eriod p - -0.03 log P + 1.39. Similarly, 
ISabbev et al.l (Il995h using hydrodynamic models to measure the 
p-factor as a function of pulsation phase, found p - 1.3 to 
p - 1.6. Recent theoretical measurements based on spectral 
line profiles from hydrodynamical spherically-symmetric model 
stellar atmospheres, with HARPS spectroscopic measurements, 
[Nardetto et al. (2004, 2007) yielded both a geometric (their Eq. 
5) and a dynamic Pp relation, assuming that radial velocities are 
determined using the first-moment method. On the other hand, 
iNardetto et al.l (l2009h measured a steeper slope for the dynamic 
relation from spectra when radial velocities are determined us- 
ing a cross-correlation technique, p - -0.08(+0.05)logP -i- 
1.31 (±0.06), hence referred to as the cross-correlated hydrody- 
namic Pp relation. The cross-correlation method for measuring 
radial velocities lead to a Pp relation that is about 5% smaller 
than the Pp relation based on the first moment method. This sug- 
gests that the p-factor is a function of wavelength and is also sen- 
sitive to the te chnique used to me asure radial velocities. It should 
be noted that Gray & StevensonI (12007 ) proposed a toy model for 
fitting the pulsation velocity to spectral line profiles directly to 
avoid the need of a p-factor, though this result is correlated to 
limb-darkening relations employed by the authors. When their 
predictions of pulsation and radial velocity are converted to a 
p-factor for the same observations as previous works, they find 
larger values of the p-factor, by about 5%. 

On the other hand, the projection factor has also been mea- 
sured by using the infrared surface brig htness techn i que to 
observe Cepheids with known distances. iGieren et al.l (l2005l) 
required the period-projection factor relation to be /? - 
-0.15 logP-H 1.58 to measure a dista nce to the LMC that i s inde- 
pendent of Cepheid pulsation period. IStorm et all (1201 Ibh found 
similar results for LMC Cepheids, also using the infrared surface 
brightness t e chniqu e, for different observations. Furthermore, 
IStorm et aP (12011 a') used observations of Galactic Cepheids 
with known HST parallaxes to verify the steeper period depen- 
dence of the projection factor; we refer to this type of analysis 
as t he LMC-HST Pp relat ion and p-factor. lMJraniiet al.1 (l2005h 
and iGroenewegenI (l2007h measured the p-factor for Galactic 
Cepheids with known distances using interferometric and spec- 
troscopic observations, both authors measured p - 1.27. One 
might argue that this method does not measure the actual value 
of the p-factor, but measures a "fudge" p-factor that corrects for 



uncertainties in the distance measurement. Another challenge is 
that angular diameter measurements imply the radius is defined 
at a layer in the photosphere that is the same as the layer that 
defines radial velocity observations; this layer may differ 

Observational and theoretical measurements of the p-factor 
disagree, leading to a significant problem if one wishes to use 
the Baade-Wesselink method to calibrate the Leavitt Law as well 
as understanding the metallicity dependence of the Leavitt Law. 
The uncertainty of the p-factor is one of the greatest sources 
of uncertainty for calibrating the Leavitt Law. The differences 
between observational and theoretical measurements are probes 
of unknown and uncertain physical effects in Cepheids. In this 
work, we consider one possible source of the differences, i.e. 
stellar limb darkening, which is one of the most important con- 
tributors for determining the p-factor Specifically, we explore 
how stellar limb darkening from hydrostatic model stellar at- 
mospheres, described in the next section, varies as a function 
of fundamental parameters, themselves defined as functions of 
pulsation period from the literature. We explore how these in- 
tensity structures depend on the assumed geometry of the model 
atmospheres in Sect. 3. In Sect. 4, we compute the period projec- 
tion factor (Pp) relation from the predicted intensity profiles and 
compare to previous results, and hypothesize why the LMC-HST 
and cross-correlated hydrodynamic Pp relations differ in Sect. 5. 
We summarize the analysis in Sect. 6. 

2. Model stellar atmospheres 

Model stellar atmospheres are computed for both plane-parallel 
and spherically-symmetric geometries. We compute s pheric ally- 
symmetric models using the S At las code (.Lester & NeilsonI 
20081). This code is based on the iKurucd d 19791) Atlas code 
but updated for geometry in Fortran 90. Plane-parallel model 
atmospheres are computed using our Fortran 90 version of the 
code. These codes assume the stellar atmosphere is in local ther- 
modynamic equilibrium and hydrostatic equilibrium. Radiation 
transfer is computed as a function of wavelength and intensity 
profi les are c omputed for 1000 /i-points for increased preci- 
sion. i Kuruczl ( 1979) mode ls pr edict intensities for 17 or less yu- 
points, while Claret (|2008|) and lClaret & BloemenI ([201 1 ) com- 
pute limb-darkening laws from intensity profiles containing 100 
//-points. The large number of yu-points are necessary to ade- 
quately compare predicti ons an d observed intensity profiles (e.g. 
Neilson & Lester 2008a[l20Tlh . 

The model stellar atmospheres are described by various fun- 
damental parameters, a plane-parallel model by specifying the 
effective temperature and gravity, while a spherically-symmetric 
model requires an additional parameter such as mass, radius, or 
luminosity. We compute model atmospheres as a functio n of pe- 
riod. The gravity of a Cepheid is given by the relation (iKovacsl 
I2OOOI) 



logg = -1.21 logP-H 2.62, 



(2) 



while the radius is given by the period-radius relation determined 
by Gieren et al. (1999), 



logiRIRo) = 0.680 logP -I- 1.146. 



(3) 



The bolome tric luminosi ty is determined by the bolometric 
Leavitt Law (lTurnedl2010h . 



log(L/L0) = 1.148 logP-H 2.415, 



(4) 



and by combining the bolometric luminosity and radius, the ef- 
fective temperature is given as a function of period. 
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We note that using these relations impHcitly assumes a 
period-mass relation. This is contentious because Cepheid 
masses determined by different met hods yield di f ferent values, 
i.e. the Cepheid mass discrepancy ( Kelleill2008t iNeilson et alJ 
1201 lah . However, this does not affect the results, as the atmo- 
spheric structure is insensitive to small (10 - 20%) changes of 
stellar mass. 

For this work, model stellar atmospheres are computed as 
a function of pulsation period for both geometries, assuming 
a range of pulsation periods from log P - 0.4 - 2 in steps of 
0.1, covering the known period range for Galactic Cepheids. For 
each model atmosphere, intensity profiles are computed for the 
BVRIH, and /T-bands, though for most of this work, discussion 
is restricted to V and /T-band profiles. It is these t wo wavebands 
that are use d for in terferometric observations (iKervella et al.l 
[2004; Moura rd et al.l B009) as well as f or the infrared surf ace 
brightness technique dGieren et al.ll2005HStorm et al.ll2011allbh . 
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3. Intensity profiles 

We explore how intensity profiles vary as a function of pulsation 
period. The V- and /T-band intensity profiles, normalized with 
respect to the stellar flux, are shown in Fig. [T]f or both geometries 
as a function of pulsation period. Intensity profiles from plane- 
parallel model stellar atmospheres do not vary much, especially 
in the V-band. Conversely, intensity profiles from spherically- 
symmetric model stellar atmospheres vary significantly because 
of the atmospheric extension of the model atmospheres that in- 
creases as a function of pulsation period. This atmospheric ex- 
tension cannot be modelled by plane-parallel atmosphere codes. 
Plane-parallel intensity profiles tend to have smaller central 
intensities than spherically-symmetric model intensity profiles. 
A model described by some pulsation period will have a spe- 
cific value of luminosity and flux that must be conserved for any 
intensity profile regardless of geometry, i.e. for the same stellar 
flux. 



I Ipv(M)f^d^i = I /sph(/u)//(i//. 



(5) 



Because spherically-symmetric model atmospheres account for 
atmospheric extension then the intensity towards the stellar limb 
for a spherical model is smaller than that for a plane-parallel 
model atmosphere, i.e. /sphC/^ — » 0) < /pp(// — » 0). Thus, to 
conserve emergent stellar flux then the central intensity for the 
given spherical model must be greater than the central inten- 
sity for plane-parallel model atmosphere described by the same 
fundamental parameters. This is equivalent to noting that in 
spherically-symmetric models the flux diverges as a function of 
optical depth, leading to a steeper temperature gradient, hence 
a sharper decrease of intensity as a function of yU. The mean in- 
tensity, J = J Iip)diJ, along with the central intensity, increases 
as well. This difference is clear in the bottom panels of Fig. [1] 
where we show the difference between V- and /T-band flux- 
normalized intensity profiles from spherically-symmetric (sph) 
and plane-parallel (pp) model atmospheres. A/ = /sph(A') - ^pp(A') 
as a function of pulsation period. The difference is greatest near 
the stellar limb, depending on the atmospheric extension of the 
spherical model atmospheres, whereas towards the center of 
the stellar disk the spherically-symmetric intensity profiles are 
brighter than the plane-parallel model atmospheres. 

We can take this a step further and define the variable 



o-i = 



Ksph >'p| 



(6) 



'pp 



Fig. 2. The relative total intensity difference between spherical 
and plane-parallel model intensity profiles as a function of pul- 
sation period, for the V-band (solid line) and AT-band (dashed 
line) intensity profiles. 

Unlike the stellar flux, the mean intensity is not required to be 
the same for stellar atmospheres of differing geometry, but de- 
scribed by the same fundamental parameters. Therefore, ctj is 
a relative measure of the model atmosphere's deviation from 
plane-parallel geometry, i.e. the amount of atmospheric exten- 
sion. We plot this fractional difference between the mean inten- 
sity for each geometry, cri, as a function of pulsation period for 
V- and /T-band intensities. At short periods, log P x 0.4, the dif- 
ference between the mean intensities for spherically-symmetric 
and plane-parallel model atmospheres is approximately 20%, re- 
flecting the amount of atmospheric extension for a stellar atmo- 
sphere with \ogg ^ 2.1. At logP - 2, the difference is almost 
50%. 

The results of Fig. |2] suggest that Cepheids cannot be 
physically represented by plane-parallel mode l stellar atmo- 
sphere s and plane-parallel radiative transfer. "Mareng o et al.l 
(20021 l2003l) computed h y drody namic stellar atmospheres us- 
ing the ISasselov & Ragal d 19921) code, radiati ve transfer wa s 
computed using the plane-parallel Atlas code (lKurucz|[l979l) . 
We showed in Fig. [T]that the intensity distribution from plane- 
parallel model stellar atmospheres do not vary significantly 
with respect to fundamental parameters, whereas spherically- 
symmetric models do. In this section, one can conclude that 
intensity profiles from plane-parallel model atmospheres are 
insensitive to assumed fundamental parameters, hence the 
pulsation period. On the other hand, intensity profiles from 
spherically-symmetric model atmospheres are more centrally- 
concentrated with increasing pulsation period, i.e. the amount of 
atmospheric extension depends on the pulsation period. This is 
important because the projection factor depends on the structure 
of intensity profiles, hence how the projection factor depends 
on pulsation period also depends on the assumed geometry of 
model atmospheres. This also suggests that limb darkening is a 
more significant effect in measuring Cepheid angular diameters 
and p-factors than found previously. 

4. Angular diameter corrections 

iMarengo et al.l (l2002h found that limb darkening varied signif- 
icantly as a function of pulsation phase, by an amount greater 
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Fig. 1. Predicted intensity profiles normalized with respect to the Eddington flux for plane-parallel model stellar atmosphere (top) 
and spherically-symmetric model atmospheres (middle) in V-band (left) and /T-band (right) wavebands. The bottom panels show 
the difference between intensity profiles computed from spherically-symmetric and plane-parallel model atmospheres. Spherical 
models predict much smaller intensities near the stellar limb relative to the plane-parallel models. For spherical models, intensity 
profiles for short-period Cepheids have smaller central intensities and sharp declines at smallest value of /i. Each line represents a 
step of A log P = 0. 1 from log P = 0.4 - 2. 



than the variation shown in Fig. [T]for differences in fundamental which assumes the intensity at any point on the star is equal 

parameters. Limb-darkening corrections, ^ud/^ld, however, do to the centra l intensity and the actual limb-darkened angular di- 

not change (Mare ngo et al. 2003) . Limb-darkening corrections ameter, 0ld (iDavis et aDl2000h . lKervella et al. (2004) presented 

are the ratio between the uniform disk angular diameter, 0ud, interferometric observations of nine Galactic Cepheids to cal- 
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Fig. 4. The predicted projection factor as a function of wave- 
length for the P = 10 day model atmosphere for the spherically- 
symmetric (solid line) and plane-parallel (dashed line) geome- 
tries. 



Fig. 3. Predicted angular diameter corrections 9\jd /6'ld as a func- 
tion of pulsation period computed from spherically-symmetric 
(solid lines) and plane-parallel (dotted lines) model stellar atmo- 
spheres. In the plot, the K-hand corrections are greater than the 
V-band corrections. 



ibrate the surface-brightness relation, however limb-darkened 
angular diameters presented depended on theoretical limb- 
darkeni ng corrections from Claret (2000) . Neither Claret (2000) 
nor Ma rengo et al.l (l2003l) considered the role of geometry in 
predicting limb-darkening relations or limb-darkening correc- 
tions. Because the central intensity of a spherically-symmetric 
model atmosphere is greater than that for a plane-parallel model 
for a given stellar flux, then 9ud for a spherical model will be 
smaller, meaning a smaller limb - darken ing correction. This sug- 
gests that iMarengo et al.l (l2002i l2003h overestimated the limb- 
darkening corrections for Cepheids as well as how much pulsa- 
tion affects the intensity profile as a function of pulsation phase. 
In Fig.[3l we present angular diameter corrections as a func- 
tion of pulsation computed from plane-parallel and spherically- 
symmetric model stellar atmospheres for V- and /T-bands. 
Corrections computed from plane-parallel model atmospheres 
do not vary by more than 0.3% for all periods, but the 
spherically- symmetric limb-darkening corrections vary from 
k = 0.93 to 0.88 in the V-band and k = 0.98 to 0.92 in 
the /T-band. These corrections vary significantly, both from 
plane-parallel corrections and as a function o f puls ation pe- 
riod. These results suggest that iKervella et al.l (l2004l) underes- 
timated the angular diameters of the observed Cepheids, and 
may explain the small but significant dif ferences between sur- 
face brightness relations measured by Kervella et al.l (l2004l) 
and Fouque & Gieren (1997). The result is also consistent with 
the re sults of ISasselov & Karovska (1 19941) and iNardetto et al.l 
(12004 . 



5. Period - projection factor relation 

Using our predicted intensity profiles, geometric p-factors (po) 
are computed using Eq. [1] as functions of pulsation period and 
wavelength. Projection factors are computed as a function of pe- 
riod as p = Pofgmdfo-g, where fg^d is a correction due to ve- 
locity gradient within the atmosphere of the Cepheid, and fo-g 
a correction due to the differential photo spheric pulsation ve- 
locity between optical and gas layers (see INardetto et al]|2007l 
for more details). We assume that /gnd - -0.02 log f + 0.99, 



and/„_^ = -0.023 log P + 0.979 from FNardetto et al.] (l2007b. We 
then multiple our Pp relation by 0.95 to correct for the depen- 
dence of the p-factor on the the cross-correlation method. This 
computation of the p-factor allows us to explore if differences in 
ste llar limb-darke ning between the models computed here and 
by INardetto et al.l (2009) may resolve the differences between 
theoretical hydrodynamic Pp relations and those derived using 
the IRSB technique. If the difference remains then we need to 
consider other uncertainties and phenomena. 

We start by considering the p-factor as a function of wave- 
length for the logP = 1 model stellar atmospheres assum- 
ing both plane-parallel and spherically-symmetric geometries, 
shown in Fig. |4] The p-factor differs significantly as a function 
of geometry where those from plane-parallel model atmosphere 
are 3-7% larger. It also varies as a function of wavelength: for 
the spherically-symmetric log P - I model stellar atmosphere 
the p-factor varies from about p - 1.22 at A x 400 nm to about 
p = 1.37 as /I — > oo. This is consistent with the differences 
between V- and A'-band intensity profiles shown in Fig. [T] At 
shorter wavelengths, the intensity profile is more centrally con- 
centrated than intensity profiles in the near infrared, suggesting 
a smaller p-factor. 

The next step is to compute the p-factor as a function of 
pulsation period. This is done for BVRIH and /T-band inten- 
sity profiles for both model geometries, and is shown in Fig. |5] 
The plane-parallel model geometric p-factors are constant with 
respect to pulsation period, variation of the p-factor seen in 
Fig. |5]is due to variation of fo-g(\ogP) and fgradi^ogP). The 
spherically-symmetric model p-factors decrease as a function 
of period non-linearly. In this case, we are comparing geomet- 
ric p-factors from our hydrostatic plane-parallel and spherically- 
symmetric model atmospheres to p-factors for four cases. The 
first case is the Pp relation from plane-parallel model atmo- 
spheres coin£^ed_to_spectra_^om Hindslev & Bell (1986), but 
also see lParsonsI (11972) and Burki et aLj (Il982l) . T he second case 
is Pp relations from spherical mod els, Eq. 6 in Nardetto et al.l 
(l2007l) . while the third case is the INardetto et al.l (l2009l) rela- 
tion. The former relation uses the first moment method, while 
the latter assumes radial velocities are measured using the 
cross-correlation method. The fourth case is the LMC-hHST 



period-projection factor relation from IStorm et aP (|20TTal!b), 
though we could also incl ude results fro mLMerand et al. (20051), 
'GroenewegenV2007') and'Ngeow et al.'(T012). The LMC+HST 
Pp relation from Storm et al. (201 la) has a slope measured from 
LMC Cepheids while the zero-point is measured from Galactic 
Cepheids. Our results appear to agree with the INardetto et al.l 
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Fig. 5. The p-factor predicted from plane-parallel and 
spherically-symmetric model atmospheres as a function of pul- 
sation period. Projection factors computed from plane-parallel 
models are shown for the BVRIH and /iT-bands as blue dashed 
lines, while projection factors from spherically-symmetric 
models are shown as red solid lines. For refe rence the period - 
proje ction factor relations a re shown from Hindslev & Belli 
(IT986) (magenta dotted line). iNardetto et al.l (l2007l), their Eq~6 
(orange dashed-double d ot line), Nardetto et al.l (l2009b (pale 
blue dot-dashed line) and IStorm et all (1201 Ibh . (black double- 
dashed line). The projecti o n facto rs for Galactic Cepheids 
determined by IStorm et aLl (1201 lab are also included (black 
dots). 



Table 1. Best-fit period-geometric projection factor relations 
from spherical model atmospheres. 
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± 0.0003 


1.414 


± 0.001 


0.0017 




H 


-0.0226 


± 0.0003 


1.438 


± 0.001 


0.0017 




K 


-0.0217 


± 0.0003 


1.444 


± 0.001 


0.0017 








Power-Law:po 


= x(\ogP) 


-+>' 








<L 




y 


z 


cr 


B 


-0.0224 


±0.0016 


1.357 


± 0.001 


1.51 ±0.08 


0.0010 


V 


-0.0236 


± 0.0014 


1.384 


± 0.001 


1.67 ± 0.07 


0.0011 


R 


-0.0239 


± 0.0014 


1.400 


± 0.001 


1.72 ± 0.07 


0.0012 


I 


-0.0247 


±0.0015 


1.418 


± 0.001 


1.74 ± 0.07 


0.0013 


H 


-0.0268 


±0.0016 


1.441 


± 0.001 


1.78 ± 0.07 


0.0014 


K 


-0.0262 


±0.0017 


1.448 


± 0.001 


1.76 ± 0.08 


0.0014 



Table 2. Parameters for Galactic Cepheids with model stellar 
atmospheres 



(|2009|) relation, but not the observed relation, verifying that 
limb-darkening does not explain the difference between the ob- 
served and theoretical relations, hence we must consider differ- 
ent p henomena to explain the difference between the relations 
from lNardetto et al.l (l2009l) and lStomi etaP (l2011al) . 

The LMCh-HST Pp relation is based on observations of 
lower-metallicity Cepheids, while our models are computed at 
solar metallicity. If stellar limb-darkening varies significantly as 
a functio n of metallicity, then the Pp relations might be resolved, 
however, Nardetto et al. (2011) determined from hydrodynamic 
models that metallicity does not resolve that difference. We com- 
puted another set of model stellar atmospheres following the 
same period-radius, period-gravity and period-bolometric lumi- 
nosity relations, but at LMC metallicity. At lower metallicity, the 
p-factor is marginally larger than at solar metallicity; by at most 
3% for our P = 100 day Cepheid. Only for the longest-period 
Cepheid does limb darkening differ significantly. We note that 
these models test only differences in metallicity and not differ- 
ences in pulsation period due to differences in metallicity. To test 
that effect, we would require a LMC period-bolometric luminos- 
ity relation (the other two relations are believed to be metallicity 
independent), which is not available. This suggests metallicity 
differences between LMC and Galactic Cepheids do not lead to 
significantly different p-factors and agrees with Nardetto et al. 
(1201 Ih . 

The predicted geometric p-factors calculated from 
spherically-symmetric models are fit with various relations. We 
assume a linear, quadratic, and power-law relations. The best-fit 
coefficients for each relation are presented in Table [T] The 
predi cted linear geometric r elations are c onsistent with re sults 
from iHindslev & Belli (IT986.) as well as iParsonsI (Il972h and 



Name 


Period 


T-eff 


R 


logg 




days 


{K) 


(Re) 


(cgs) 


^Cep 


5.37 


5956 


42.6 


2.0 


j/Aql 


7.18 


5772 


49.0 


1.8 


XCyg 


16.39 


5283 


100.6 


1.15 


SVVul 


45.00 


5224 


168.7 


0.7 



i Nardetto et al.l (|2007|) (their Eq. 5 for the geometric p-factor). 
The quadratic relations do not fit the predicted projection factors 
much better than a linear model, but the power-law fits best 
with an exponent ranging from 1.5 to 1.75. Thus, geometric 
p-factors, hence p-factors themselves are non-linear functions 
of period. 

As a consistency check, we compute model stellar atmo- 
spheres for a sample of Galactic Cepheids, listed in Tab. |2] with 
measured gravities, effective te mperatures and radii , derived us- 
ing the IRSB technique from IStorm et aP (1201 lah . Projection 
factors for these Cepheids are computed as a function of wave- 
length, which we compare to our predicted Pp relations in Fig.|6] 
The predicted p-factors for each Cepheid are similar to the pre- 
dicted Pp relations for each waveband, however the Cepheid p- 
factor is sensitive to the measured gravity of the Cepheid. For in- 
stance, if we compute a model with the same parameters as SV 
Vul but with logg = 1.15 then the p-factor increases by about 
3-5%. The sensitivity of the p-factor on the assumed gravity also 
grows with decreasing gravity, i.e. increasing pulsation period. 
This may help explain some of the difference between the LMC- 
HST and cross-correlated hydrodynamic Pp relations. 
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Fig. 6. Projection factors predicted from spherically-symmetric 
model atmospheres representing a sample of Galactic Cepheids 
as a function of waveband, from B to A'-band. In the plot, the 
colour green represents B-band, blue is V, magenta is R, light 
blue is /, orange is H and black is A'-band. 



6. Why do theory and observations disagree? 

There is a clear difference between LMC-HST and cross- 
correlated hydrodynamic Pp relations, with no obvious resolu- 
tion. We explore three possibilities: Cepheid mass loss and in- 
frared excess; the effect of Cepheid chromospheres; and possible 
dispersion of the Pp relation due to the possible range of Cepheid 
gravities. The difference between the two types of Pp relations 
may also be related to the uncertainty o f radial veloci t y mea sure- 
ments due to line profile asymmetries. 'Ngeo w et al.l (1201 2h dis- 
cusses how differences in radial velocity measurements change 
the determined p-factor for 6 Cephei. In terms of the Pp relation, 
these differences will not lead to a systematic shift of the Pp re- 
lation but will lead to a greater scatter of the p-factors used to fit 
a relation. 

Recent observation s s how that Cepheids undergo mass loss 
(Marenao et al. 2010ab; B armbv et aP l20lU iMatthews et al.l 
120 11; .Neilson et al...201 lb) . This mass loss is believed to lead 
to a far-infrared ex cess and may b e consi s tent with the K-band 
exces s observed by 'Kervella et alj (l2006l) : iMerand et aP (l2006i 
120071) . [Neilson et al. (2010) showed that LMC Cepheids ap- 
pear to undergo mass loss based on me asured optical and IR 
fluxes compi l ed by .Ngeow et al.l (l2009h . From this analysis, 
iNeilson et al.l (1201 Oh showed that Cepheid mass loss affects the 
observed {V - K)q colors, hence the infrared surface brightness 
method is also affected. In that case, the IRSB technique predicts 
a 3 - 4% larger angular diameter than reality, at least for shorter 
periods. This leads to the Pp relation that is also 3 - 4% overes- 
timated. This difference is not enough to explain the difference 
between cross-correlated hydrodynamic and LM C-HST Pp rela- 
tions, but more detailed observations are needed (iKervella et al.l 
1201 Ih . 

Pulsation-driven chromospheric activity is another possi- 
bility. Phase-dependent chromospheric activity le ads to limb- 
brigh tening near the edge of the stellar disk (Sasse lov & Lesten 
Il994a,b.) . Consider a Cepheid with a simple limb-darkening re- 
lation 



from iGettiiigl (11934 . which leads to a p-factor of p = 1.41. If 
we assume limb-brightening near the edge of the stellar disk 61 
then the geometric p-factor becomes 



where we assume that terms J dlj^^dii x; 0. As a result. 



S9{ 

Pnew = 1-41(1 + ---r- 

ri{) 



(8) 



(9) 



m , 3 

= 1 _ (1-^), 

Kn =1) 5 



(7) 



If p„,„ = 1.5 then d^HMa = 0.064. This suggests a 6% 
increase of stellar flux due to a change of intensity near the 
limb of the stellar disk. Assuming a limb-brightening relation 
6I(jS)/I{fi = 1) = fl, where a is constant when /u < 0.05 and 
zero elsewhere, then for a = 20, S'H/'Hq - 0.06. It is un- 
likely that a Cepheid chromosphere can produce the necessary 
intensities to resolve this difference in the p-factor, but may still 
contribute. Pulsation-generated chromospheres are believed to 
be generated by shocks and shock pr operties depend on pulsa - 
tion period and effective temperature (iNeilson & Lesterll2008bl) . 
At short period, shocks generated by pulsation are stronger and 
more frequent than at longer period, furthermore. X-ray and UV 
observations of short-period Cepheids suggest the presence of 
very hot > 10^ ^ plasma in the envelope that may be related to 
chromospheres (Engle et al. 2009). This simple test suggests that 
chromospheric effe cts may affect the Pp relation, for instance, 
ISabbev et al.l (1 19951) did find a value p - 1.6 from hydrodynamic 
models suggesting the possibility of limb-brightening in spectral 
lines as a function of pulsation phase. 

A third possible resolution is related to uncertainties of 
Cepheid parameters, in particular gravity. We showed in the pre- 
vious section that a change of A log g = 0.35 for S V Vul changes 
the predicted p-factor by about 3-5% depending on wavelength. 
Cepheids can have a large range of gravity based on the range 
of known radius from about 30 - 200 Rq, but also because of 
their masses. Ceph eid masses ar e uncertain; this is the Cepheid 
mass discrepancy ('Keller '2008^ and may also be affected by 
Cepheid mass loss (Neilson et al. 201 la). The combined disper- 
sion of masses and radii as a function of pulsation period leads to 
a signifi cant dispersion i n Cep heid limb darkening, hence the Pp 
relation. iNardetto et al.l (|2007|) measured the Pp relation from a 
sample of eight Cepheid models. The predicted relation may be 
affected by the small sample that does not account for the disper- 
sion of the Pp relation due to differences in gravity for the same 
pulsation period. 

Unfortunately, none of these three phenomena are obviously 
sufficient to uniquely resolve the differences between cross- 
correlated hydrodynamic and LMC-hHST Pp relations, but a 
combination of the three effects might do so. This analysis also 
highlights just how sensitive a precision understanding of the p- 
factor is to the physics describing Cepheids. 



7. Summary 

The Cepheid projection factor is one of the most significant 
sources of uncertainty for calibrating the Cepheid Leavitt law 
to the precision necessary to constrain cosmological parameters. 
This uncertainty is compounded by the fact that various meth- 
ods for determining the projection factor and its dependence on 
pulsation period do not agree. This uncertainty is also affected 
by errors in previously published angular diameter corrections 
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used for interferometric observations. Plane-parallel model at- 
mospheres predict limb-darkening corrections that are 3 to 6% 
larger than spherically-symmetric model atmospheres. This dif- 
ference will be explored in greater detail in a future article. 

In this work, we computed p-factors from intensity pro- 
files predicted using spherically-symmetric model stellar atmo- 
spheres as a function of wavelength. Each model is computed 
for assumed values of luminosity, radius and gravity, each a 
function of pulsation period. From the models, the Pp relation 
is predicted for a number of wavebands and compared to the 
results from numerous works, but do not agree with observed 
relations (IStorm et al.ll2011 a). However, the predicted projec- 
tion factors are consistent with measured projection factors for 
Galactic Cepheids with astrometric distances. 

We also explored three possible reasons why there is 
such significant disagreement between cross-correlated hydro- 
dynamic and LMC-HST Pp relations. Cepheid mass loss and 
infrared excess may affect the observed colors, in turn adding 
uncertainty to the predicted angular diameters measured using 
the IRSB technique, or cross-correlated hydrodynamic models 
do not account for potential limb-brightening due to pulsation- 
driven chromospheres. Conversely, both LMC-HST and cross- 
correlated hydrodynamic relations are derived by different sam- 
ples, where the measured p-factors are very sensitive to the grav- 
ity of Cepheids, which is in turn, has a significant dispersion as 
a function of pulsation period. None of these three hypotheses 
appear to resolve the differences, but could contribute to a por- 
tion of the differences. This does highlight the sensitivity of the 
p-factor to atmospheric dynamics and Cepheid properties. 

In summary, we conclude that spherically-symmetric model 
stellar atmospheres provide insight into understanding the pro- 
jection factor, the Pp relation and the wavelength dependence 
and dispersion of these relations. These models also hint at a 
possible nonlinear Pp relation that reflects the nonlinear nature 
of Cepheid atmospheric extension as a function of period. This 
non-linearity may be more pronounced if we employed a non- 
linear bolom etric Leavitt Law, s imilar to the nonlinear optical 
Leavitt Law (iKanbur et al.ll2004l) . 
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